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Abstract 

Technical aspects of the Shirkov-Solovtsov's analytic perturbation the- 
ory (APT) are considered. We construct explicitly two sets of specific 
functions, {2t„(s)} and {A n (Q 2 )} that determine the nonpower asymp- 
totic expansions for Minkowskian and Euclidean QCD observables in 
APT. The results, up to third order, are written in terms of the Lam- 
bert W-f unctions. As an input we used the exact two loop and the three 
loop (corresponding to Pade transformed beta-function) RG solutions for 
' common invariant coupling a s . In addition, the exact three-loop cou- 

pling is expanded in powers of the exact two- loop solution. The excellent 
accuracy is achieved with few terms of this series. We derive order by or- 



(N 



der elegant systems of equations for both sets of the functions. Then we 



o 
o 



43 



construct the global versions of the APT functions with quark thresholds 
in the M S scheme and give numerical results. 



1 Introduction 

In recent publications |2j a modified version of QCD perturbation theory, free of unphysical 
t^j- 1 singularities, was elaborated. It was termed the renormalization group (RG) invariant analytic 
approach ||. In this works the idea to combine renormalization invariance and Q 2 -analyticity 
for the running coupling was explored. New important properties of the analytic coupling 
established. This makes it possible to gain some information about infrared behavior of the 
Oh! theory starting from perturbative calculations. 

In this paper we shall be concerned with analytic perturbation theory (APT) suggested 
in Ref. ||], a particular version of RG invariant analytic approach. A thorough mathematical 
investigation of the nonpower APT series has been undertaken in works ||, where the stabil- 
ity of APT results (as compared to the conventional perturbative approach) with respect to 
the renormalization scheme and higher loop dependence for the whole low energy region was 
explained (see also ||). APT has been successfully used for description of many important 
spacelike and timelike processes (for a partial list of contributors to this subject, see ref. H). 

For the reader's convenience here we briefly summarize some aspects of APT. Let D(Q 2 ) 
be a dimensionless quantity, depending on a single energy scale Q. A suitable example for our 
purposes is Adler D-function related to some timelike process. Usually, it is presented in the 
form of power series f] 

oo 

D P t{Q 2 ) = A)(i + d n<{Q\ /)), (i) 

n=l 

here Dq is a process dependent constant and / denotes the number of active quark flavors at 
the energy scale Q. In APT, the same quantity should be presented in the form of a nonpower 
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asymptotic expansion as 

oo 

D an (Q 2 ) = A>(i + d nMQ 2 J)) (2) 

n=l 

where A n is the "analyticized n-th power" of the coupling in the spacelike region: A n (Q 2 , f) = 
{a™(Q 2 , f)} an - Let R(s) be the physical quantity determined through D{Q 2 ) in the timelike 
domain (for example R e + e -(s) or R T (M T )). Then in APT it has the representation the || |7j 

R(s)=R (l + r{s)) r(s) = J2Zidn%n(s,f). (3) 

The functions 2l n are defined through the transformation 

1 f s+le dz 

%n(sj) = -— -A n (-Z,f), (4) 

2m J s _ le z 



previously introduced in Refs. || [| (see also recent paper [[UJ). The inversion of ([4[) reads 

r°° d? 

A n (Q 2 J)=Q 2 J o {s + Q2)2 *n(s,f). (5) 

In recent works [||, [7] the universal functions A n and 2l n have been calculated analytically at the 
one loop order. In the infrared region, an oscillating behavior for these functions was observed. 
To construct these functions, beyond the one loop order, the iterative approximation for the 
running coupling was used |7|, [11], [12| or RG equation was solved numerically in the complex 
domain 1131. 



Recently exact two-loop order solution to the RG equation has become available [14], [15 
The solution has been written in terms of the Lambert-W function. In addition, the third 
order RG equation (with Pade transformed /3-function) has also been solved in terms of the 
same function ji~5|| . These explicit solutions proved to be very useful for determination of the 
analyticity structure of the coupling in the complex Q 2 -plane. Using the exact two loop solution 



the analytically improved coupling Af\Q 2 , f) was reconstructed |L4|, |16| . 

Afterwards, in paper |T7||, a higher order solution to the RG equation (in arbitrary MS-like 
renormalization scheme) was expanded in powers of the scheme independent explicit two-loop 
order solution. In this way, a new method for reducing the scheme ambiguity for the QCD 
observables has been proposed. A similar expansion for the observable (depending on a single 



energy scale Q) motivated differently, has been suggested in yl 

The aim of this paper is to give practical formulas for calculation of the APT sets {A n (Q 2 , /)} 
and {2l n (s, /)}. As an input we use the above mentioned explicit solutions to the RG equa- 
tion. In Sec. 2 we derive order by order general equations for both sets of the APT functions. 
In Sec. 3 we present exact two-loop results. The Minkowskian functions are constructed, in 
the closed form, in terms of the Lambert-W function. The corresponding spacelike functions, 
A n (Q 2 , f), are reconstructed through the spectral representation, or equivalently using inverse 
transformation (|5|). In Sec. 4 the APT functions are calculated to the three loop level. The 
Pade approximation improved three-loop coupling is used. General method for constructing 
the APT sets, beyond two-loop level, is discussed in Sec. 5. In Sec. 6 we construct global 
universal functions, A n {Q 2 ) and 2l n (s), (both introduced in 0). For crossing the quark flavor 
thresholds, we use the continuous matching conditions for the MS scheme coupling a s (Q 2 , f). 
Sec. 7 contains the numerical results. Some technical details are given in the Appendix. The 
conclusions are given in Sec. 8. 
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2 General results 

Let us settle conventions and notations. The running coupling of QCD satisfies the RG equation 
Q 2 das ^ 2 ,f) =P f MQ 2 J)) = -jtplap>{C?,f), (6) 

n=0 

with the normalization condition a s (p 2 , f) = g 2 /(47r), where p is the renormalization point. In 
the class of schemes where the beta-function is mass independent (3q and (3{ are universal 



the result for /3| in the modified MS" (MS) scheme is |]TP 

, _ 1 /2857 5033/ 325/ 2 \ 

^-(i^V~2 ir + ^r)- (8) 

The Euclidean functions are defined through the spectral representation 

A n (Q\ f) = - T ^Qda = - r , t 9/ ^ Pn(t, f)dt, (9) 



vr J v + Q 2 vr J _^ (e* + Q 2 /A 2 

where the spectral function p n (a, f) = ^s{a s (— o — «0)} n , A is the QCD scale parameter, and 
p n (a, f) = p n (t, /) with t = ln(cr/A 2 ). The timelike set of functions {2l„(s, /)} is defined by 
elegant formula 

1 f°° dcr 

% l (sJ) = - —Pn(<T,f), (10) 

Js cr 



obtained in work [E0 



In numerical calculations singular integral @ should be regulated. With the help of (0), 
integral (Q) can be represented as 

MQ 2 , f) = MQ 2 , f, T) + 2l n (A 2 e T , /) + 0(e" T ), (11) 

where A n (Q 2 ,f,T) denotes the integral (0) taken over the finite interval —T < t < T. For 
sufficiently large values of T (when exp(— T)(Q 2 /A 2 ) <C 1), the contributions of order e~ T can 
be neglected. So that the second term, 2l n (A 2 e T , /), compensates the main truncation effects in 
the integral. Formula ( |TTD enables us to achieve a good numerical precision even for moderate 
values of the cutoff T. 

In what follows, we shall derive the equations for the functions A n (Q 2 ,f), 2t„(s, /) and 
p n (a, f). In the k-loop order, the results have the form 

a f£' /) = -"Efl^+'W'-fl- n=1.2..., (12) 
^^ = -nJ2p f N ^l N+1 (sJ), n = l,2.... (13) 
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^C^ = -E^U% +1 (*,/), n = l,2.... (14) 

The similar one-loop order equations were obtained in paper 0] . We start from the RG equation, 
in the k-loop level, presented as 

^^o^-^^r^'w 2 ,/) (15) 

where, n > 1, is some integer. The APT variant of Eq. ([T^) reads 

\ da l^oP ) =-«E^W 2 >/). (16) 

^ ^ J an at — n 



N=0 



on the other hand, we can write 



da n s {Q\f)\ 1 r d ° ^daU-a-tOJ) 

3 KT-7 _ ~7v\ ' ( 17 ) 



aing 2 j an 7t7 a + g 2 ain(-<7-io) 

From the identity, ln(— a — zO) = In a — Z7r, it follows that 



3- 



da^-a-zOJ) = d%a?(-cr-iO,f) = dp n (cr, f) 
<91n(— a — iQ) din a din a 

so that 

da^Q 2 ,f)\ 1 [°° da dp n (a,f) 



dlnQ 2 ) an irj a + Q 2 dlna 



(19) 



an 

integrating by parts fll9|) and using the asymptotic vanishing of p n (cr, /) (the asymptotic free- 
dom) 

-^ Pn (a,f)\™ = 0, (20) 

we get 

da«(Q 2 ,f)\ 19 p rfa , f , dA n (Q 2 J) 



dlnQ* J an ndlnQ2j a + Q^ Pn ^ ^ «91nQ 2 ' (21) 

from this and (0) follows system of equations (12). From C[T~0f) , we see that system (|14D is a 
consequence of system ([[3]). Let us multiply Eq. ([L4]) by the factor (a + Q 2 ) -1 and take the 
integral over the region < a < oo. Integrating by parts and taking into account condition 
(p0|), we recover the system (|T2|). Note that for n = 1 Eqs. ( [T2]) and ( |T3| ) are analogous to the 
basic Eq. @ with a" replaced by 2t„ and .A n respectively. We remark, that Eqs. flUMED are 
derived on general grounds using the RG equation (^) and the spectral representation together 
with the asymptotic freedom condition (f20|). 



3 Exact two-loop results 



For convenience, in what follows we shall omit index f in the coefficients f3^. The exact two-loop 
solution to Eq. (|) is given by [[14], [T^ 



here b\ = Pi/Pq, A = and W(C) is the Lambert W function [[Hp, the multivalued inverse 
of 

C = ^(C)exp^(C). 

The branches of W are denoted Wk(C), k = 0, ±1, .... 

By performing analytical continuation of function ( p2] ) and its powers in the complex Q 2 - 

plane f], we determine the corresponding spectral functions pn\a, f) = f>n\t, /), n — 1, 2 

For < / < 6 0, the result is 



Pi 2) (t,/)=(|)^ 



1 + Wi(z(t)) 



with 



*(*) = 7— exp(-t/6i +i(l/6i - 1)tt). 
Oie 



(23) 



(24) 



After insertion of (p3|) into (|Tl] ) and (0) we construct explicitly the two-loop functions .An ^ (Q 2 , /) 
and 2l£ (s, /). The integrals for the Minkowskian functions can be performed \ We obtain 



^ 2) (^/) = ^lnfl + |«i 2 )(- S )) 



24 2 W) 



.CJ^l n (l + ^ a (2)(_ a ) > \_ §L a W( 



A) 



ft 



(25) 
(26) 
(27) 



/) = (-£)" ±* {in (l + |a»(-*)) - |«-(- S ) + , e tcm 

where, for convenience, we introduced the notation 

with z s = ^(s/A 2 ) -1 ' 6 ^ 1 ^ 1 / 61-1 ) -1 . Note that, the functions 2l n (s, /) are determined through 
(n — 2)-th order residual terms of the Taylor expansion of ln(l + j-a^(-s)) in powers of 

ft« (2) (-*)- 

Using the asymptotic properties of the W-function [BT|], in the limit s — > 0, we recover the 
result of work [^D| 

2tS 2) ( S ,/)^^, and 2#W)^0 /or n > 1. (30) 

Note that, the "analyticized" nonpower perturbative expansions for the timelike observables, 
may be rewritten as power series in traditional coupling a s (s), but with modified by 7r 2 -factors 

d For details of analytical continuation we recommend papers ji~4| , and . 
e For / > 6 formula (||) should be changed, see @, [jlj. 

f The one-loop expressions for the Minkowskian functions were derived in Refs. |t], |2^| . 
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coefficients. Let us consider the two loop case. The function ( pop has the following formal 
power expansion 



where a s (s, f) is the exact two loop coupling (22). Higher order functions 2l n 



n 



s, /)+..., 
(31) 
1,2...) have 



similar power series expansions. Substituting these expansions into (|3|) we find the expansion 
for .R(s) in powers of a s with modified by 7r 2 -factors coefficients. Previously, in the timelike 
region, similar expansions in powers of the iterative solution to RG equation, were obtained 
in IIHI]. Application of the series can be found in papers |23|, |24|, 25]. The "7r 2 -effects" for 
various timelike quantities have been estimated in paper [25|, in particular, it was found that 
the 7r 2 -factors give dominating contributions to the coefficients of R e + e -{s) = a to t(e + e _ — > 
hadrons)/<r(e + e _ — > fi + p~). On the other hand, in recent papers p6| , p7[ various timelike 
events were analyzed in the f=5 region. Higher-order "7r 2 -effects" have been taken into account 
properly. It was found that the extracted values for a s (M 2 ) are influenced significantly by these 
effects. 



4 The three-loop case with Pade approximation 

Pade improved perturbative series attracted much interest recently. They, as opposed to trun- 
cated perturbative series, exhibit reduced renormalization scale and scheme dependence ||28|| . 
Pade related resummation method which achieves elimination of the unphysical dependence 
from observables was proposed in recent works [p9| . 

The solution to RG Eq. at the three-loop order, with Pade transformed beta-function 

(3 Pade = -Poa't I 1 + ' " 



Po — pT a s 



has the form 15 



(3) ,rii t\ A) 1 ,1 „.._ {0A\ ( Q 2N _fr 



a ^ (QJ) = -Ja-M^ + w. l( iY i = -eb-r p {w)\v) ■ < 32 > 

Starting from (|52j), one can construct Apl din (Q 2 , f), just like in the two-loop case. For < 
/ < 6, the spectral function, p Pa de,n{cr, f) = p Pa de,n{t, f) = 3{a^(-cr - i0)}, is defined as 

^ Jt,/) = (l) *\T^ww+wm) ' (33) 

with 

Z(t) = — expiPofo/ff - t/bx + <l/6i - 1)tt). (34) 
Integral ([10|) , with the spectral function (|3~3"D , can be done. Thus we obtain 

4 3 L,i(^ /) = "4- (^HWi(Zs)) + (1 - -)QHV + Wi(Z s )) - tt) , (35) 
^Po \V V J 
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*f^M) - ± ~ d " b* ) ) . (3a) 



where 77 = 1 — P0P2/P1 and 

Z s = -(s/A 2 )- 1 ^ exp(-r ? + .(l/^ - l)vr). 


For n > 2, the following recurrence formula is valid 

(37) 

where p n is given in (|55|). Putting together (|3~6|) and (|37D we get 



3^ _ ^ 



( n _l)( r/+W / 1 (^))n-l 

5 Multi-loop case 



n > 2. (38) 



As follows from the paper |T7| one can obtain the expressions for multi-loop functions in terms 
of the two loop functions. In particular, the coupling of higher orders can be represented as 

00 

^ fe) =E c W n ' ( 39 ) 

n=l 

with Ci = 1. The analyticized (Euclidean and Minkowskian) versions of fl3~9|) read 

00 

4 fe) = £ c ^ 2) > (40) 



n=l 



^ = J2 ci nV- (41) 
n=l 

So the two-loop functions can be considered as the minimal basis of any orders perturbative 
expansions ]17|]. Any observable O (except the quantity possessing an anomalous dimension) 



can be represented as the series 

00 

OW = J>«,e (42) 

n=l 

We would like also to emphasize that the one-loop coupling function could not be used for a 
similar expansion: the case is that multi-loop coupling functions have more complicated singu- 
larities f| However the exact two-loop coupling (expressed through the Lambert W-function) 

s For example, the two-loop solution contains double logarithms In In a; that evidently can't be expanded in 
powers of 1/ Ins. 
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can be used to reflect correctly the higher orders contributions. Note that some similar results 
was obtained in work |18| . 



By using the limiting conditions (|30| ) in the representation (|4T|), we derive immediately 
universal infrared limit |L[ for the Minkowskian coupling to any finite order 

»?W)-4-, s^O. (43) 

Formula (|43"1) has been obtained in work |20| in a different way. Evidently, the same limit has 
the Euclidean coupling A^\Q 2 } f) at Q 2 — > (the remarkable result of works So we can 

express any observable in any loop order in terms of the Lambert-W function. The expressions 
obtained will have correct analytic properties and finite infrared limit. 



6 The quark thresholds 

In MS-like renormalization schemes important item is how to introduce the matching conditions 
at the heavy quarks thresholds for the strong coupling constant . For pedagogical introduction 
to this subject we recommend paper |3D| (see also Ref. [3~I 1). In the context of APT, this 
problem was studied in works J7|, |26], |27j. The matching conditions and analyticity have been 
combined there via the special model spectral function 

Pn(a) = p{= 3 (a, A 3 ) + M 2 )(p{(a,A f ) - p{ t ~ 1 (a,Af^ 1 )), (44) 

/>4 

where the mass Mf corresponds to the quark with flavor / and the values for A f are determined 
through the continuity conditions 

a s (M 2 ,f) = a s (M 2 ,f-l) f = 4,5,6. (45) 

To be more exact, in the MS-like renormalization schemes (beyond the one loop level) conditions 



(fiq) should be modified |3(| • With the modified conditions the final results are not sensitive to 
the exact scale one uses to connect the couplings. Nevertheless, in this work we assume simple 
formula ( f45| ) up to third order: the case is that in APT this simplification does not introduce 
numerically significant errors. Note that with the explicit solutions (^) and (|32|) , the equation 
(|45|) for Af can be solved analytically. For instance, inserting the three- loop solution (^) in 
([4~5|) we solve 

A 7 = Mf{-b f F(z f ^) expfa/ + F(z f ^))} b f/ 2 (46) 



where b f = (3[/(^) 2 , Vf = 1 - / (p{) 2 



J f-i \ M f 



A f -N 2/ ^ 



Pt* Pi 



F{zf-i) = (Vf-i + W.^Zf.,))^^ - Vf . (4, 



Using ( |4^ ) with (|46| ) we construct global quantities A n (Q 2 ) and 2l n (s) relevant in the whole 
ranges of Q 2 and s respectively. 
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7 Numerical results 



The numerical calculations are performed by using the computer algebra system Maple V 
release 5; Maple has an arbitrary precision implementation of all branches of the Lambert- W 
function. 

(3) 

In Table 1, various 3-loop approximants are compared with exact three- loop coupling, atnum, 
the numerical solution of the transcendental equation 

, Q 2 1 &i , f 1 h \ 2&a - b\ ( b x + 2b 2 (3 a s b x \ 

In —77 = — In ^» - + — h b 2 H f= — arctan = arctan — = , 

A 2 A,« s 2 V/3o«l A>« fl V VA V VA y/Aj 

(49) 

where = fik/ (3q +1 (k = 1, 2) and A = 4b 2 — b\. By we denote the commonly used three 
loop order iterative approximant [31] 



" ' ' = /%! ~ /f ~P~ U (1 ° ' ftj' (50) 

where L = InQ 2 / A|j . Usually the same formula fl5Tj) is used in the timelike region. a ts 
denotes the truncated series (|3~9"D 

^-E^, (51) 



in this case cf } = 1, cj, 3) = 0, 4 3) = cf } = and cf ] = §/3 2 //? 2 . From the Table 1, we 

see that, a ts allows us to achieve more good agreement with the exact answer than the Pade 
and iterative approximants. a ts is accurate, for Q > 900 MeV (in these region it practically 
coincides with the exact answer) while the Pade approximant worsens below Q = lGOOMeV^ 
and iterative approximant allows us to achieve 1% accuracy only for Q > 2300 MeV. We have 
taken A 3 = 400 MeV, and the values of Af (/ = 4, 5, 6) are determined by imposing conditions 
(53) on the exact numerical solution fl4"9p. For the quarks masses throughout this paper we 
assume the values M x = M 2 = M 3 = 0, M 4 = 1.3 GeV, M 5 = 4.3 GeV and M 6 = 170 GeV. 

In Table 2 the results for the Minkowskian expansion functions are given in the three 
quark region (0.4 GeV< <2.6 GeV). The functions 2l^ de>n (s, /), %£\sj) and Sl£J n (s, /) 
are compared for n = 1,2 and / = 3. We take the common value, A/ = 3 = 400 MeV. The 
differences between and 21^ are less than 2.19%, while the differences between 2l[ 3 \ and 
%-Pade i are smaner than 0.12%. For the second functions (n = 2), these differences are of the 
same order: \8% {2 \%)\ < 2.6 and | 8%f ade2 (%)| < 0.3. 

Table 3 summarizes the results for Euclidean expansion functions, An\Q 2 , /), in the two 
loop and three-loop orders, for n — 1, 2 and / = 3. The three quark region is chosen (0.4 GeV 
< Q < 2 GeV) and we take Af = $ = 400 MeV. The numerical calculations have been performed 
using formula ( jTT| ) (high precision was achieved already for T = 1000). In the three loop case, 



(3) (3) 

two approximants, A\ sn and A PadSn , are compared. There is good agreement between these 
approximants ( 0.07% end better for n=l, and 0.14% and better for n = 2). The differences 
between A\ sn and A n (n = 1, 2) are larger, up to 1.72% for n = 1 and up to 1.4% for n = 2. 

The asymmetry fSJ, 5 as (%) = (A%\(Q 2 , f) - Slg^Q 2 , f))/A£] x (Q 2 , f) * 100, was found to 
be appreciable (compare Tables 2 and 3): it increases from 2% at Q = y/s = 0.4 GeV to 7.5% 
at Q = v/i = 2.0 GeV. 
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On the other hand, for global expansion functions A n {Q 2 ) and 2l n (s) (that are constructed 
according the model flSID ) all the above mentioned differences are slightly enhanced. In Table 5, 
the results obtained with Minkowskian global approximant, W ts l (s), are given. For comparison 
we include the differences, S 2 \%) = (Slg^s) -2lf ) (s))/2tg > ) 1 (sj * 100 and d (3 \%) = (21^ (s) - 

% ( pade,i( s )) /2£i(s) * 100. As an input we take A /=3 = 400 MeV. The values of A f (/ = 
4, 5, 6) for various approximants (determined through (flop) are given in Table 4. The maximal 
differences are observed at « 1.3 GeV: S 2 \%) = 3.21 and S 3 \%) = -0.63 respectively. 
The analogical results for the global Euclidean coupling Ai(Q 2 ) are given in Table 6. In this 
case DW(%) < 2.64 and |D( 3 )(%)| < 0.5 for 0.4 GeV < Q < 100 GeV. 



8 Conclusions 

In this paper we have constructed analytically two sets of the specific functions, {A n } and 
{2l n }, that determine the APT expansions for the spacelike and timelike observables respec- 
tively. The results, up to third order, are written in terms of the Lambert W-function. For the 
both sets of the functions, we have derived order by order the infinite sets of differential equa- 
tions. These equations follow from the RG equation, spectral representation and asymptotic 
freedom. With the help of the equations, at the two-loop level, compact expressions for the 
Minkowskian functions 2t n (s, /) (see Eqs. (p5"D - (p8|) ) are obtained. For the Euclidean functions, 
A n (Q 2 , f), the regulated expressions (|TT| ) guaranteeing high precision in numerical calculations 
are presented. We have shown the connection between APT series for the timelike observables 
and the corresponding power series with modified (by "7r 2 -factors" ) coefficients (see Eq. (EH])). 
For Pade- improved 3-loop case, the Minkowskian set of functions {2t n (s, /)} is constructed 



see fl38Q ). In addition, to calculate the APT functions in the three-loop order, the general 



method of work |T7| is used. The advantage of this method over the method of Pade approxi- 
mation is demonstrated (see Table 1). 

The matching conditions for the running coupling function a s (Q 2 , /) at the flavor thresholds 
to three loops in MS scheme are presented in analytic form (see (|46"D ). Global functions A n (Q 2 ) 
and 2l n (s) with flavor thresholds are constructed. 

We have examined numerically some of the APT functions to second and third orders (see 
Tables 1-6). The accuracy of various approximants to the spacelike and timelike APT functions 
is estimated. 

Let us discuss the status of the obtained exact and approximate expressions for the functions 
A n and 2l n . The question may be raised as to whether they allow us to achieve advantage 
over the existing in the literature iterative formulas. In recent works [|K], [11, |i~2f , analytical 



continuation to the timelike region has been performed using the iterative solution (|50| ) to the 
RG equation. However, it was estimated in ref. J39|, that with the iterative solution, in the 
two loop case, the error in A^(Q 2 , /) and %l[ (s, /) may be as large as 4%-5%. Furthermore, 
the iterative formulas are not convenient from the technical point of view too. Since, starting 



from iterative formulas, we have not found analogous of Eqs.(12)-(13) which greatly simplify 
numerical analysis. On the other hand, the precision of the experimental data (for R e + e - (s), R T , 
ets.) is continuously increasing. Therefore the authors hope that the more accurate theoretical 
formulas will be helpful for practical calculations in the APT framework. It should be remarked, 



that an alternative to our method was given in work [[L3|l , where the RG equation, to third order, 
has been solved numerically in the complex domain. 

Other possible applications of the obtained formulas are in the context of different (not 
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APT) approaches to Landau pole problem |33|, 
frameworks developed in works |36| [37], [38|] . 
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34, 35] and more general (non-perturbative) 



9 Appendix 



In this appendix we calculate analytically 2l n (s, /) to second order. Let us define the auxiliary 
functions 

1 



Rn(sJ) 



7T 



dta n (tj). 



(52) 



111 s 



where s = s/A 2 and Qa n (t,f) = p n (t, f) = p n (cr, /) with a = exp(i). Then 2l n (s, /) 
QR n (s, f). The expressions for a, can be read from (|23f) . In the two loop case 



l 



/3 1 i + w 1 (z(t)y 

With ( |53"1) integral (]52[) can be rewritten as a contour integral in the complex z-plane 

dz 1 



(53) 



R { n\sJ) 



Pn 



z + 



(54) 



where, we denote 



Pr. 



vr/3i V A 



ra-2 



2 S = z(lns), 



z f = ee u i ; , 



(55) 



and the limit e — > is assumed. Let us change the integration variable in (|54|), = W(z), we 
then get 

—— —duj. (56) 

■f lW ^(l + c^) n 



For n > 2, from (1561) we obtain the equation 



01 



(n-2) 



n = 3,4... 



(57) 



the two-loop version of the equation fllcf) . It is sufficient, to calculate 21^ and Sli, 2 '*. Indeed, 
higher functions can be determined using recurrence formula fl5"7|). Note that R{ (s, f) is 
divergent in the limit e — ^ 0, (see fl54|) ). Nevertheless, it has finite imaginary part []. By direct 
calculation we obtain formulas (g5) and (Effi), formulas ( p7|) and fl28|) follow from formula ( ^q) . 

h For the asymptotic behaviour of the W function see paper [Ell. 
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Table 1: The various approximations to the three loop coupling versus exact numerical solution 
otnlm of the RG equation. We include the relative errors 5a[f(%) = (ai« m — a[f)/anum * 100, 
£«p 3 L(%) = (ail - af ade )/a^ um * 100 and 8af '(%) = (a£L - a^/a&L * 100. 



Q GeV 


a£L(Q 2 ,f) 


8af s \%) 




8af(%) 




/ = 3 


A 3 = 400 MeV 




.8 


.764911 


1.90 


-18.88 


-15.5 


.9 


.633231 


.82 


-7.42 


-9.3 


1.0 


.554140 


.43 


-4.40 


-6.1 


1.1 


.500275 


.26 


-3.04 


-4.3 


1.2 


.460747 


.17 


-2.28 


-3.3 


1.3 


.430253 


.12 


-1.81 


-2.6 


1.4 


.405868 


.09 


-1.49 


-2.2 


1.5 


.385829 


.06 


-1.26 


-1.9 


1.6 


.369009 


.05 


-1.09 


-1.6 


1.7 


.354645 


.04 


-.96 


-1.5 


1.8 


.342204 


.03 


-.85 


-1.3 


1.9 


.331300 


.03 


-.77 


-1.2 


2.0 


.321646 


.02 


-.70 


-1.1 


2.2 


.305269 


.02 


-.59 


-1.0 


2.4 


.291842 


.01 


-.51 


-.9 


2.6 . 


.280585 . 


.01 . 


-.45 


-.9 




/ = 4 


A 4 = 354.407 MeV 




2 


.330496 


.009 


-.40 


-.39 


4 


.244183 


.002 


-.15 


-.29 


6 


212689 


.001 


-.10 


-.26 


8 


195076 


.0004 


-.08 


-.24 


10 


183390 


.0003 


-.06 


-.23 




/ = 5 


A 5 = 259.602 MeV 




10 


.187687 


.00003 


-.02 


-.01 


20 


.160349 


.00001 


-.01 


-.04 


30 


.147853 


.00001 


-.01 


-.04 
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Table 2: Comparison of the Minkowskian functions 2l„ (s, /), %i Pade n (s, f) and 2Q,/ n (s, /)■ We 
denote = (2lg? n - 2l^ 2) )/2lg^ * 100 and 0ag>* >B (%) = (<n " 41e,J/<n * 100. 







<52l 2) 


(%) 


°^Pade,l 


(%) 2li 3 Us 3) 




u<K Pade,2\ /Q ) 








f — 
J — 


3 A, 


— 4DD MpV 






4 




1 97 




1 1 
. 1 1 


. louooj 


1 QQ 


98 


c; 




1 fi^ 




D7 


1 994fi^ 


i 

-l.UO 


.ou 


.u 


49^Q7 


1 QD 




.uo 


1 1490^ 


1 8 


9Q 


7 


.oyooz i 


z.uo 




ni 

-.Ul 


. luuuiy 


.Ol 


9^ 


Q 

.o 


^771 1 n 

.O i i 11U 


9 14 




D4 




1 04 


91 


.9 


.359032 


2.18 




-.06 


.093863 


1.45 


.16 


1.0 


.343733 


2.19 




-.08 


.088593 


1.75 


.12 


1.2 


.319204 


2.16 




-.10 


.079832 


2.16 


.05 


1.4 


.300337 


2.11 




-.11 


.072913 


2.39 


-.01 


1.6 


.285311 


2.04 




-.12 


.067343 


2.51 


-.04 


1.8 


.273011 


1.98 




-.12 


.062776 


2.57 


-.07 


2.0 


.262720 


1.91 




-.12 


.058965 


2.60 


-.09 


2.2 


.253955 


1.85 




-.12 


.055737 


2.60 


-.10 


2.4 


.246378 


1.79 




-.12 


.052967 


2.59 


-.11 


2.6 


.239747 


1.74 




-.11 


.050561 


2.57 


-.12 



Table 3: Comparison of the Euclidean functions, Ah (Q 2 ,f), at the two-loop and three loop 
orders. 



Q GeV 


< j (Q 2 ,3) 






4 2, (Q 2 ,3) 


•ApideaiQ 2 ' 3) 


^ts^iQ 2 ' 3) 


.4 


.507853 


.511591 


.511933 


.118913 


.117136 


.117273 


.6 


.438444 


.443587 


.443761 


.107177 


.106285 


.106436 


.8 


.393408 


.399067 


.399119 


.097030 


.096780 


.096911 


1.0 


.361380 


.367164 


.367133 


.088705 


.088884 


.088990 


1.2 


.337219 


.342947 


.342862 


.081875 


.082339 


.082421 


1.4 


.318218 


.323809 


.323686 


.076208 


.076860 


.076922 


1.6 


.302807 


.308224 


.308076 


.071443 


.072219 


.072264 


1.8 


.290003 


.295233 


.295069 


.067385 


.068241 


.068273 


2.0 


.279159 


.284203 


.284028 


.063886 


.064794 


.064815 


2.2 


.269831 


.274693 


.274511 


.060839 


.061779 


.061790 


2.4 


.261702 


.266390 


.266205 


.058160 


.059117 


.059121 


2.6 


.254538 


.259063 


.258876 


.055784 


.056749 


.056747 
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Table 4: The flavor dependence of the MS scheme scale A/, as determined from the approxi- 
mants a4um, a[f and a® rf -. 



f 


3 


4 


5 


6 


A^ 2) MeV 


400 


341.444 


242.195 


102.113 


Ail,/ MeV 


400 


354.407 


259.602 


112.194 


A® MeV 


400 


354.034 


259.288 


112.043 




400 


358.667 


263.932 


114.306 



Table 5: Comparison of various approximants for global Minkowskian coupling 5li(s). d^(%) 
and d^(%) denote relative errors for 21^ (s) and 5lp^ respectively. 







d^{%) 


d®(%) 




y/s GeV 


o2i(«) 


d^{%) 


dW(%) 


.4 


.520162 


1.98 


-.23 




3 


.245069 


2.61 


-.54 


.5 


.477242 


2.41 


-.30 




4 


.224780 


2.42 


-.51 


.6 


.443950 


2.70 


-.36 




5 


.211545 


2.27 


-.48 


.7 


.417373 


2.90 


-.42 




6 


.201948 


2.17 


-.46 


.8 


.395663 


3.03 


-.47 




7 


.194467 


2.08 


-.45 


.9 


.377585 


3.11 


-.51 




8 


.188409 


2.01 


-.43 


1.0 


.362286 


3.16 


-.55 




9 


.183363 


1.95 


-.42 


1.1 


.349158 


3.19 


-.58 




10 


.179068 


1.90 


-.40 


1.2 


.337757 


3.21 


-.60 




20 


.155103 


1.63 


-.36 


1.3 


.327751 


3.21 


-.63 




25 


.148685 


1.55 


-.34 


1.4 


.318641 


3.17 


-.62 




30 


.143821 


1.50 


-.33 


1.6 


.303211 


3.08 


-.61 




50 


.131743 


1.36 


-.30 


1.8 


.290593 


2.99 


-.60 




70 


.124838 


1.28 


-.28 


2.0 


.280039 


2.91 


-.59 




90 


.120137 


1.23 


-.27 


2.4 


.263277 


2.77 


-.57 




100 


.118270 


1.20 


-.27 
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Table 6: Comparison of various approximants to global Euclidean function Ai(Q 2 ). The relative 
errors D^(%) and D®(%) correspond to A?(Q 2 ) and A { ^ Pade (Q 2 ) respectively. 



Q GeV 






DW(%) 




Q GeV 








.4 


.530379 


1.50 


-.25 




7 


.205366 


2.22 


-.46 


.6 


.462088 


1.97 


-.32 




8 


.198339 


2.15 


-.45 


.8 


.417301 


2.25 


-.37 




9 


.192520 


2.10 


-.44 


1 


.385158 


2.42 


-.41 




10 


.187592 


2.05 


-.43 


1.2 


.360725 


2.53 


-.44 




20 


.160570 


1.75 


-.38 


1.4 


.341394 


2.59 


-.46 




30 


.148141 


1.60 


-.35 


1.6 


.325635 


2.62 


-.48 




40 


.140458 


1.50 


-.33 


1.8 


.312489 


2.64 


-.49 




50 


.135046 


1.44 


-.32 


2 


.301321 


2.64 


-.49 




60 


.130940 


1.39 


-.31 


3 


.263297 


2.58 


-.50 




70 


.127669 


1.35 


-.30 


4 


.240696 


2.48 


-.50 




80 


.124975 


1.32 


-.29 


5 


.225362 


2.38 


-.48 




90 


.122699 


1.29 


-.29 


6 


.214096 


2.29 


-.47 




100 


.120740 


1.26 


-.28 
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